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Late time decay of very heavy dark matter is considered as one of the possible explanations for 
diffuse PeV neutrinos observed in IceCube. We consider implications of multimessenger constraints, 
and show that proposed models are marginally consistent with the diffuse 7-ray background data. 

Critical tests are possible by a detailed analysis and identification of the sub-TeV isotropic diffuse 
7-ray data observed by Fermi and future observations of sub-PeV 7 rays by observatories like HAWC 
or Tibet AS-I-MD. In addition, with several-year observations by next-generation telescopes such as 
IceCube-Gen2, muon neutrino searches for nearby dark matter halos such as the Virgo cluster should 
allow us to rule out or support the dark matter models, independently of 7-ray and anisotropy tests. 


PACS numbers: 95.35.-l-d 95.85.Ry, 98.70.Vc 

The origin of cosmic high-energy neutrinos [in is a 
new mystery in astroparticle physics (see, e.g.. Refs. [3- 
[^). Various theoretical interpretations include possibili¬ 
ties of hadronic {pp) production in cosmic-ray (CR) reser¬ 
voirs 0 and photohadronic (py) production in hidden 
CR accelerators [IMl , and the observed neutrino in¬ 
tensity at ~ 0.1-1 PeV energies is consistent with earlier 
models [iB - ll^ . Only a fraction of the observed events 
could have Galactic origins (e.g.. Refs, (lol-l^). 

Not only astrophysical sources but also dark matter 
may lead to high-energy neutrinos and 7 rays (see recent 
reviews, e.g., Refs. [^, [^). Because of several motiva¬ 
tions such as the thermal relic hypothesis and unitarity 
bounds [23 - [^ , most studies had focused on dark matter 
with TOdm ^ 30-100 TeV. However, there is no funda¬ 
mental objection to considering very heavy dark matter 
(VHDM), which is hard to probe by existing accelerators 
such as the Large Hadron Collider. As considered prior 
to the IceCube observation, indirect searches in neutrinos 
and 7 ra ys give us unique opportunities to high-energy 
searches . Assuming nondetections of cosmic neu¬ 

trino signals, in light of IceCube and Fermi, the power 
of multimessenger approaches had been demonstrated to 
constrain particle pro perties of VHDM [29l - l^ , even for 
Wdm ^0.1 PeV [H, ll^. As soon as PeV neutrinos 
were discovered, the VHDM scenario was invoked [ssl - 
[S^ and various phenomenological models have been de¬ 
veloped [ssl - ld^ . Although they do not give a natural ex¬ 
planation why the observed neutrino flux is comparable 
to both the diffuse 7 -ra y ba ckground and CR nucleon- or 
nuclei-survival bounds , the VHDM scenario can 

presently be consistent with the data [H, . 

In order to test various possibilities, the multimessen¬ 
ger approach and point source search are essential. Their 
power has been demonstrated in Refs. [§,[H,li3illIl and 
Refs. [52l - [^ . respectively. In this work, we consider 
how these two strategies can be used to test the VHDM 


scenario with current and future observations. 

The VHDM Scenario. — The mean diffuse neutrino 
(and anti-neutrino) intensity is calculated by evaluating 
line-of-sight integrals. Although we calculate it numeri¬ 
cally throughout this work, for decaying VHDM, the all 
flavor intensity is analytically estimated to be 
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where and are extragalactic and Galactic con¬ 
tributions to the cumulative neutrino background, re¬ 
spectively (e.g.. Ref. [13 )■ The VHDM decay scenario 
predicts similar Galactic and extragalactic contributions. 
We have used h Ri 0.7, Utm ~ 0.3, Da « 0.7, ildmh^ = 
0.12, pcC^ = 1.05x 10“®/i^ GeV cm“^, tn is the age of the 
Universe, Psc<? = 0.3 GeV cm“^ in the Solar neighbor¬ 
hood, and Rsc = 8.5 kpc. Note that ~ 0.6 corrects for 
redshift evolution of decaying VHDM [H, 113, and Jq, 
is the dimensionless J factor averaged over D [ 2 ^, [13 . 
We use the Navarro-Frenk-White profile to show results, 
but for decaying VHDM we checked that our basic con¬ 
clusions are not altered for more cored profiles. Pre¬ 
dictions for the diffuse 7 -ray intensity and single source 
fluxes should be very similar, since their normalization is 
fixed by the diffuse neutrino intensity. The VHDM life¬ 
time Tdm = 7dm,27.5 10 ^^'^ s is a model parameter to be 
constrained, and TZv = Tl,z{E,j) is the energy-dependent 
function converting the bolometric flux to the differential 
flux at E,,, which depends on final states (e.g.. Ref. [13 )■ 
Assuming that all decay products are Standard Model 
particles, for demonstration, we consider several models 
proposed by Refs. [13[130- Following Refs, [s^ [53- 
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FIG. 1: Diffuse all-flavor neutrino and 7 -ray intensities ex¬ 
pected in the VHDM scenario. The ES13 model is assumed 
with Tdm = 3.0 X 10^^ s. The total (thick dashed line) and 
extragalactic (thin dashed line) contributions to the cumula¬ 
tive neutrino background are shown with the observed data. 
The expected 7 -ray background is also shown (thick solid) 
with the latest Fermi data. We also show contributions of 
extragalactic cascaded 7 rays and direct 7 rays from Galac¬ 
tic VHDM, which are not affected by uncertainty of Galactic 
magnetic fields. KASGADE and CASA-MIA 7 -ray limits are 
indicated. 

with electroweak corrections, the final state spectra ob¬ 
tained from 10 TeV to 100 TeV masses are extrapolated 
to PeV masses. Our choice of VHDM models is such that 
they include both hard and soft spectra, so our results 
can be viewed as reasonably model independent [ 2 ^ [ 2 ^ . 

In Figs. 1 and 2, we show examples of the viable 
VHDM scenario for diffuse PeV neutrinos observed in 
IceCube. Using the ES13 model [s^, where the VHDM 
mass TOdm = 3.2 PeV is used, we consider DM —s- 
and DM ^ qq with 12% and 88 % branching fractions, 
respectively. Although a bit larger masses are favored to 
explain the 2 PeV event, one can easily choose param¬ 
eters accounting for the observed data. In the RKP14 
model [4l|, the Majorana mass term is introduced in the 
Lagrangian, which may lead to metastable VHDM de¬ 
caying into a neutrino and Higgs boson. Reference 
suggested another interesting scenario, where the light¬ 
est right-handed neutrinos constitute dark matter with 
TOdm = 0{1) PeV. We also consider this model for 
TOdm = 2.4 PeV, assuming branching fractions DM ^ 
: DM ^ vZ : DM —>■ « 2 : 1 : 1, where the 

neutrino spectral shape turns out to be similar to that of 
Ref. [4l| (see Fig. 2). As in the latter two models, spec¬ 
tra may be more prominently peaked at some energy, and 
VHDM does not have to explain all the data. 

7-Ray Limits. — Standard Model final states from 
decaying or annihilating VHDM lead to 7 rays as well as 
neutrinos. If final states involve quarks, gluons and Higgs 
bosons, neutrinos largely come from mesons formed via 
hadronization, and 7 rays are produced. A spectral bump 
is produced by two-body final states such as uh and/or 
weak bosons via leptonic decay into a neutrino and 



FIG. 2: The same as Fig. 1, but for the RKP14 model with 
Tdm = 3.5 X 10^^ s. 

charged lepton. Electroweak bremsstrahlung is relevant 
even for possible decay into neutrino pairs. In extragalac¬ 
tic cases, the fact that the diffuse neutrino and 7 -ray in¬ 
tensities are comparable gives us generic limits 
In Galactic cases, 7 rays below ^ 0.3 PeV can reach the 
Earth without significant attenuation, air-shower arrays 
such as KASGADE and GASA-MIA as well as 
Fermi [ 6 l| provide us with interesting constraints [l^ . 

We numerically calculate the diffuse 7 -ray background, 
including both extragalactic and Galactic components. 
Thanks to the electron-positron pair creation, suffi¬ 
ciently high-energy 7 rays are attenuated by the extra¬ 
galactic background light and cosmic microwave back¬ 
ground. Then, the pairs regenerate 7 rays via the inverse- 
Compton and synchrotron emission. For an extragalac¬ 
tic component, we calculate electromagnetic cascades by 
solving Boltzmann equations. The resulting spectrum 
is known to be near-universal, following a C omp tonized 
E~‘^ power-law in the 0.03-100 GeV range [s^. For a 
Galactic component, it is straightforward to calculate 
primary 7 rays that directly come from VHDM. The 7 - 
ray attenuation is approximately included by assuming 
the typical distance of i?sc, which gives reasonable re¬ 
sults [ 13 . Extragalactic cascaded 7 rays (including at¬ 
tenuated and cascade components) and Galactic primary 
7 rays with attenuation unavoidably contribute to the 
diffuse 7 -ray background (see Eigs. 1 and 2). In addition, 
electrons and positrons from VHDM [9^ make secondary 
7 rays via inverse-Compton and synchrotron emission in 
the Galactic halo, as included in Eigs. 1 and 2 assuming 
a magnetic field strength of 1 /iG. Our results would be 
conservative, and weaker magnetic fields can somewhat 
increase 7 -ray fluxes. Eor cascade components, the re¬ 
sults are not sensitive to detailed spectra of final states 
from VHDM decay. See Ref. for technical details. 

Clearly, 7 -ray constraints are powerful. In the sub- 
PeV range, while the VHDM models are still allowed, 
the expected diffuse 7 -ray intensity can slightly violate 
the existing sub-PeV 7 -ray limits from old CR-induced 
air-shower experiments such as KASGADE. Thus, as we 
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here show, the proposed VHDM models can be criti¬ 
cally tested by near-future TeV-PeV 7 -ray observations 
with the High-Altitude Walter Cherenkov Observatory 
(HAWC), Tibet AS-I-MD, and perhaps by Fermi. Our 
results show that the Galactic direct component should 
be dominant above TeV energies. The VHDM scenario 
predicts that the diffuse 7 -ray intensity and large scale 
anisotropy due to Galactic components should increase 
at > 1 TeV up to 0.3 PeV, which can be tested. To eval¬ 
uate anisotropic 7 -ray emission, we calculate the factor 
averaged over the Galactic center region within 25 deg, 
and obtain Jq, ~ 6 . 8 . The excess due to Galactic VHDM 
(i.e., diffuse 7 -ray emission after isotropic emission is sub¬ 
tracted) is ~ 8 X 10“® GeV cm-^ s"! sr"! 

from Eq. (1). For emission from the Galactic plane, 
HAWG can reach ~ 5 x 10“® GeV cm“^ s“^ sr“^ 
at ~ 10 TeV and Tibet AS-I-MD will achieve ~ 
10“® GeV cm“^ s“^ sr“^ at ~ 100 TeV in five years 641. 
Hence, anisotropic TeV-PeV 7 rays from VHDM can be 
seen at least for three models considered in this work. 
Also, as clearly seen in Fig. 3, many of the diffuse neu¬ 
trinos are found in the Southern Hemisphere, outside 
the KASGADE field of view. Although diffuse TeV- 
PeV 7 -ray limits for the Galactic halo will be powerful 
enough, having 7 -ray detectors in the Southern Hemi¬ 
sphere should be much more helpful fl^ . 


In the sub-TeV energy range, extragalactic cascaded 7 
rays are relevant, and the expected diffuse 7 -ray intensity 
is marginally consistent with the Fermi data. Decompos¬ 
ing the diffuse isotropic background, although it is model 
dependent, leads to tighter constraints (cf. Refs. 
Following Ref. [^, we calculate cascaded 7 -ray bounds 
on VHDM lifetimes, using the latest Fermi data [ 6 l[. 
When total contributions are considered, we obtain lower 
limits (95% credible), = 2.3 x 10^^ s in the ES13 
model and = 1.6 x 10^^ s in the RKP14 model, 
respectively. More conservatively, for direct and extra- 
galactic contributions, we get = 1.3 x 10^^ s in the 
ES13 model and = 0.8 x 10^' s in the RKP14 model, 
respectively. The diffuse isotropic background shown in 
Figs. 1 and 2 is obtained by subtracting resolved point 
sources, so it does not involve uncertain subtraction of 
unresolved sources. Also, more than half of the total 
isotropic background in the sub-TeV range is attributed 
to resolved blazars [ 6 l|, HB, ■ Thus, the fact that the 
consistency with the 7 -ray backgrounds is marginal is 
quite robust, leading to profound implications. First, 
the diffuse 7 -ray data representing the sum of unresolved 
sources could be improved in future by Fermi, or pos¬ 
sibly HAWG. If more blazars are resolved and they give 
~ 100 % of the present diffuse isotropic background, there 
will be little room for the VHDM scenario. Second, we 
use the high-energy IceCube data presented in Ref. [B], 
which give the high significance. The extended analyses 
suggest softer spectra with a higher intensity of ^ 

10"^ GeV cm-2 s"! sr'i at ~ 20-60 TeV [Udi, but 
the lower-energy data suffer from more systematics due 
to the atmospheric muon background and possible con¬ 


tamination by Galactic sources. If they are established 
as the nearly isotropic signal, this strong case requires 
shorter lifetimes of Tdm 10 ^^ s and the diffuse 7 -ray 
background would be violated without subtraction of un¬ 
resolved sources, hinting at a different component for 
<0.1 PeV neutrinos [zB]. Our result strengthens the im¬ 
portance of understanding the < 60 TeV neutrino data. 

Muon Neutrino Limits from Galaxies and 
Galaxy Clusters. — The more direct and important 
test can be carried out by muon neutrino searches for 
nearby sources. A search for astrophysical emission is 
presented by Ref. 0 . In the VHDM scenario, the cu¬ 
mulative neutrino and 7 -ray backgrounds are dominated 
by low-mass dark matter halos, but nearby massive ha¬ 
los associated with nearby galaxies and galaxy clusters 
can be detected as point or extended sources. Following 
Ref. [IBl, we examine five nearby clusters (Virgo, For¬ 
nax, Perseus, Goma, Ophiuchus) using parameters pro¬ 
vided in Ref. [tB]. In addition, nearby galaxies M31 [tB], 
Large Magellanic Cloud (LMC) (tB] and Small Magel¬ 
lanic Cloud [zBl are considered. The signal is stronger 
for objects with larger where Mdm is the dark 

matter halo mass and d is the distance. Virgo, Fornax, 
M31, and LMC are of particular interest, and they have 
Aldm/d^ a few x 10^® Mq Mpc“^. The IceCube obser¬ 
vatory has its highest sensitivity for point source emis¬ 
sion in the Northern Hemisphere utilizing up-going muon 
neutrino events. For this reason we focus our following 
discussion on Virgo and M31. Note, however, that a pro- 
posed km® scale neutrino telescope like KM3Net in 
the Mediterranean Sea should be helpful for neutrino ob¬ 
servations from Fornax and LMC in the Southern Hemi¬ 
sphere. Although we numerically evaluate signal fluxes, 
for example, the muon neutrino flux is estimated to be 



X 


1 MdmC^ 

127rd2 TdmTlu 

1.3 X 10”®° GeV cm-2 s“® t“^_ 27.5 

/ Mdn. \ ^ 

V 5 X 10^^ Mq / V 16 Mpc / 


( 7 ^./ 15 )-^ 

( 2 ) 


for the Virgo cluster. Then, following Refs. IB, 0, 
we calculate detection rates of through-going muon 
tracks within a maximal angular range A^max — 
max[A0res, 0.5°(F,y/TeV) ' ], where the angular reso¬ 

lution is set to A0res = 0.5° for IceCube and the 
second term is due to the intrinsic uncertainty from the 
kinematics of the interaction. Although the astrophys¬ 
ical background [IBl is accounted for, the atmospheric 
backgrounds (that are taken from Refs. [tB - IsB ] ') are more 
relevant in our case. If a source is extended, can be 
regarded as the flux integrated over the source extension. 
But the backgrounds also increase, so optimization maxi¬ 
mizing the signal-to-background ratio is possible [IB, HIl ■ 
Since dark matter substructures do not play a relevant 
role for the decay scenario, the simple point source search 
is reasonable. Our results are conservative since the lim¬ 
its can be improved by analyzing starting muon tracks 
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FIG. 3: Sky distribution of showers (circles) and tracks (di¬ 
amonds) with time-ordered event numbers, with angular un¬ 
certainties. Events 28 and 32, which are likely backgrounds, 
are removed. The KASCADE field of view is shown by shaded 
regions, and a large part of the Southern Hemisphere is not 
covered. Dashed curves indicate the regions, in which 25% 
and 50% of neutrino emission from VHDM is included. Stars 
indicate positions of some nearby sources. 

and/or neutrino-induced showers for extended sources. 

In Fig. 4, we show forecasted limits that can be 
placed by searches for muon neutrinos from Virgo and 
M31. For simplicity, we assume that a next-generation 
IceCube-Gen2 detector has an effective point-source sen¬ 
sitivity that is about 5 times better than IceCube, due 
to the combination effect of enhanced effective area and 
event reconstruction [s^. We assume that this detector 
would be fully operational after the deployment season 
2019/2020, i.e., ten years after IceCube has reached its 
full fiducial volume, although quantitative results might 
be affected by details of the detector configuration. The 
90% C.L. limits are obtained based on Ref. [s^. Note 
that, although stacking analyses for nearby sources could 
improve limits in principle, we find that including objects 
with Mdm/d^ 10 ^^ Mq Mpc“^ does not help in our 
case. Their individual neutrino fluxes are too low, mak¬ 
ing the overall signal-to-background ratio worse. One 
sees the present IceCube is not large enough to test the 
VHDM scenario requiring Tdm (3-6) x 10^^ s, even with 
twenty years of operations. We need a better angular 
resolution, with which we can put crucial constraints in 
several years. This conclusion will hold for cored profiles 
even if the J factor is reduced by a factor of 2. Nonde¬ 
tections will rule out the VHDM scenario independently 
of the other limits, while positive detections may be sup¬ 
portive or suggest other astrophysical scenarios Q . 

Summary and Discussion. — The discovery of cos¬ 
mic neutrinos opens up a new window to probe new 
physics beyond the Standard Model, such as neutrino 
self-interactions [t^, and Lorentz-invariance viola¬ 

tion [sgl - l^ . The VHDM scenario has been considered as 
an explanation for the cosmic neutrinos. We considered 
two critical tests that are feasible with current and near¬ 
future 7 -ray detectors and next-generation neutrino tele¬ 
scopes. (1) The proposed VHDM models predict the dif- 



Observation Time with IceCube [yr] 

FIG. 4: Muon neutrino limits on the VHDM scenario, ex¬ 
pected for the Virgo cluster and M31. We consider the ES13 
model (solid line), RKP14 model (dotted line), and HKS14 
model (dot-dashed line), and VHDM lifetimes explaining the 
cumulative neutrino background are indicated by the shaded 
region. We assume through-going muon tracks seen in Ice- 
Gube (thick line) and a next-generation detector like IceCube- 
Gen2 (thin line) with a relative improvement of the sensitivity 
by a factor of 5. The VHDM scenario can be ruled out or sup¬ 
ported in three to five years. 

fuse 7 -ray background that is compatible with the Fermi 
data. The marginal consistency implies that they can be 
ruled out or supported by improving the data, decom¬ 
posing the sub-TeV background, and finding anisotropy 
increasing as energy. Note that the latest results of the 
IceCube Collaboration indicate a softer neutrino spec¬ 
trum with the higher intensity in the ~ 30 TeV energy 
range , which would increase the tension with 7 - 

ray bounds. (2) The diffuse sub-PeV 7 -ray background 
is also marginally consistent with the current limits. The 
excess emission around the Galactic center can be de¬ 
tected by 7 -ray and CR detectors such as HAWC, Tibet 
AS-I-MD and IceTop. (3) If the VHDM scenario is cor¬ 
rect, muon neutrinos from nearby galaxies and galaxy 
clusters such as Virgo should be detected with a next- 
generation detector such as IceCube-Gen2. Remarkably, 
this method enables us to test various VHDM models 
that only explain the data in the PeV range. 

The tests proposed here are complementary to the 
large-scale anisotropy of the arrival distribution of neutri¬ 
nos. So far, no significant anisotropy has been observed. 
We stress that our approaches become especially impor¬ 
tant if the excess around the Galactic center exists. 

Although we focused on decaying VHDM, applica¬ 
tions to annihilating VHDM are possible. The unitar- 
ity bound, which usually gives stringent limits on mdm, 
could be alleviated if the signal largely comes from sub¬ 
structures with low velocity dispersion (3^ . Although 
the predicted arrival distribution is different, constraints 
from the diffuse 7 -ray background can similarly be pow¬ 
erful. With large boost factors, muon neutrino searches 
for nearby sources are relevant as well (^ . 
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